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Abstract 



J> ' Using the Very Large Array (VLA) in the A-configuration, we have obtained a high-resolution 3.6 cm 

map of a hard X-ray source detected by the Chandra X-ray Observatory in a protostellar clump in Orion 
' molecular cloud 3. Two radio continuum sources were detected in the vicinity of the X-ray source, both 

of which have NIR counterparts. We conclude that these VLA sources are free-free emission produced by 
■ shocks in protostellar jets from the NIR class I protostars. Using the centimeter data, we determined the 

I power and orientation of the protostellar jets. The center position of the X-ray emission was found to be 

. ~ 1-2" offset from the exciting sources of the jets, and the displacement is in the direction of the jets and 

molecular outflows. We discuss the nature of the X-ray emission as the shock-excited plasma at the shock 
front where the jet propagates through interstellar medium at a speed of 1000 km s~^. 
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c/2 . 1. Introduction ular clouds 2 and 3 (OMC-2/3; Z)=450pc) with the 

^ ' Advanced CCD Imaging Spectrometer (ACIS) onboard 

Since the discovery of X-ray emission from star-forming Chandra (Tsuboi et al. 2001; Tsujimoto et al. 2002a) 

regions (SFRs) using the Einstein Observatory in the in conjunction with deep near-infrared (NIR) follow-up 

^ , 1980's (Feigelson, DeCampli 1981; Montmerle et al. 1983), imaging using the University of Hawaii 88 inch telescope 

■ low-mass young stellar objects (YSOs) have been known (Tsujimoto et al. 2003). Among 385 X-ray detections, we 

to be strong X-ray emitters. Subsequent X-ray observa- found that ~ 72% sources have a NIR counterpart while 

tories further revealed that almost all low-mass YSOs in the remaining sources have none. Based on the X-ray 

the class I — III evolutionally stages emit X-rays and the spectral and temporal properties as well as the spatial 

emission originates from optically-thin thermal plasma distribution and the K-hand luminosity function of these 

with fast rises and slow decays in flux. These lead sources, we concluded that most of the NIR-idcntified and 

to the general consensus that X-rays from low-mass NIR-unidentified X-ray sources are YSOs and background 

YSOs are attributable to high temperature (r=10^- active galactic nuclei, respectively. 

10^ K) plasma maintained by occasional flares that are Among the NIR-unidentified X-ray sources, however, 

triggered by magnetic reconnections (see a review by a few coincide spatially with the string of protostellar 

Feigelson, Montmerle 1999). clumps in OMC-3 (MMS 1-10) seen in 1.3 mm continuum 

Recent long-exposure observations by the Chandra X- emission (Chini et al. 1997). The hard X-ray emission 

ray Observatory are beginning to reveal new X-ray- (source No. 8 in Tsuboi et al. 2001; we hereafter call this 

emitting phenomena in SFRs with its unprecedented sen- source TKH 8) detected at the protostellar clump MMS 2 

sitivity and spatial resolution. These include the X-ray is one of these peculiar X-ray sources. Its associations with 

detection from a Hcrbig-Haro object (Pravdo ct al. 2001) the protostellar clump, with molecular outflows in CO, 

and the discovery of diffuse X-ray emission from high-mass IICO+ (Aso et al. 2000; Williams ct al. 2003), and shock- 

SFRs (e.g., Townsley et al. 2003). excited II2 (Tsujimoto ct al. 2002b) emission lines, and 

1.1 1 ■■ r-. ■ 1 with protostellar jets in centimeter continuum (Rcipurth 

We conducted a ~100ks observation on Orion molec- , , 1 ;i ■ i-i i . 1 11 1 

et al. 1999), make this source umikely to be a background 
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source. It is also unlikely that this X-ray emission is from 
the magnetic activity of a low-mass class I-III objects be- 
cause the X-ray emission is centered at a position that is 
significantly (~ 1-2") offset from the closest NIR sources 
(IRS 3 and IRS 5; Tsujimoto et al. 2002b). These results 
lead Tsuboi et al. (2001) and Tsujimoto et al. (2002b) 
to speculate that this X-ray emission is from a hidden 
class protostar, which is invisible by definition at NIR 
wavelengths (Barsony 1994). 

We have carried out the highest-resolution centime- 
ter imaging observation yet of TKH 8 using the Very 
Large Array (VLA) to investigate further the nature of 
this peculiar source. Given the fact that most proto- 
stars show free-free emission originating from protostel- 
lar jets (Rodriguez et al. 1995; Andre 1996), long-baseline 
interferometer imaging at centimeter wavelengths is the 
most accurate method of determining the position of jet- 
exciting sources, and the strength and orientation of the 
jet at its root. A D-configuration observation has al- 
ready detected a 3.6 cm source at TKH 8 (Reipurth et 
al. 1999), although its coarse spatial resolution is insuf- 
ficient to match the Chandra image. In this paper, we 
report the result of our A-configuration VLA observation 
of TKH 8, based on which we present another interpreta- 
tion of the X-ray emission. 
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Fig. 1. The 3.6 cm intensity map (contours) with the posi- 
tion of the NIR sources (crosses). The position accuracy of 
the VLA sources are ~0"02, while that of the NIR sources 
are ~0."1. Contour levels are 3—9<t with a step of 3cr, where 
the background noise is ~ 7.7/iJy beam~^ . The synthesized 
beam size is at the bottom left. The position angle of VLA la 
is shown with the arrow. 



2. Observation &: Result 

We observed TKH 8 with the VLA of the National 
Radio Astronomy Observatory (NRAO) on February 11, 
2002. We used the A-configuration to achieve sufficient 
spatial resolution to match that of the Chandra image. 
A 3.6 cm map was obtained with an integration time of 
~3.5 hours and the phase center at R. A.=05'^35™18?3, 
Decl.=-05°00'33" (J2000.0). The map is sensitive to 
structure smaller than ^ 2" with an angular resolution of 
~0"1. This enables us to have a much finer view of the 
X-ray emission and its vicinity compared to the prior D- 
configuration observation, which had a resolution of ~ 8" 
(Reipurth et al. 1999). 3C48 (3.25 Jy) and 0541-056 
(0.98 Jy) were used as the flux and phase calibrators, re- 
spectively. 

Data reduction, calibration and analysis were per- 
formed using the Astronomical Image Processing System 
(AIPS). The naturally- weighted map is shown in figure 1. 
We detected two sources (VLA la and VLA lb) above the 
3(7 level, for which we derived the position and the flux 
density (table 1). These two sources were not resolved 
in the prior D-conflguration observation and were named 
altogether as VLA 1 (Reipurth et al. 1999). We detected 
~ 70% of the total flux measured in the D conflguration, 
indicating that we are missing some extended emission as- 
sociated with VLA 1. VLA la is slightly extended, so we 
also determined the length of the major (minor) axis and 
position angle to be 0"40 (0"09) and 32.2 degree, respec- 
tively. The direction of the elongation is shown with the 
arrow in figure 1. 



3. Discussion 

3.1. The Origin of the Centimeter Emission 

VLA la and VLA lb are very close to the NIR sources 
IRS 5 and IRS 3, respectively (Tsujimoto et al. 2002b). 
The apparent separation between IRS 5 and VLA la 
(-O.'l) and between IRS 3 and VLA lb (~a.'07) are 
smaller than the NIR position tmcertainty of ~ 0."1. Thus, 
the NIR and radio sources are most likely coincident. Both 
IRS 3 and IRS 5 can be classified as class I protostars from 
their NIR excess and high extinction of Ay >50 mag us- 
ing the {H-K) / {K-L') color-color diagram (figure 3b in 
Tsujimoto et al. 2002b). The excess emission was also 
confirmed at the mid-infrared bands, securing the class I 
nature of these sources (Nielbock et al. 2003). 

For the following four reasons, we conclude that both 
VLA la and VLA lb are free-free emission from Hn re- 
gions ionized by the UV radiation from the shock front 
produced by a protostellar jet propagating through ambi- 
ent matter (Curiel et al. 1987). 

First, these centimeter sources are accompanied by 
class I protostars (IRS 3 and IRS 5). Embedded pro- 
tostars are frequently associated with centimeter emis- 
sion. Detailed studies indicate that most of these ra- 
dio sources, if not all, are jet-induced free-free emission 
(Anglada 1996). 

Second, IRS 3 and IRS 5 have the NIR magnitudes and 
colors of Jo > 11.3 mag and J~H > 1.6 mag, and Jq > 
11.3 mag and J-H > 3.8 mag, respectively (Tsujimoto 
et al. 2002b). This indicates that both sources have a 
mass less than 2Mq, which rules out the possibility that 
the centimeter emission is from Hn regions generated by 
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source R.A. Dccl. flux* NIR 

ID (J2000.0) (J2000.0) (/iJy) counterpart^ 

VLAla 05'^35"^I8?335 -05°00'32."97 128±18 IRS5 
VLAlb 05h35"U8f271 -05°00^33.^^93 54±11 IRS 3 

* I'he tiux density is corrected tor the primary beam response, 
t The nomenclature follows Tsujimoto et al. (2002b). 



stellar UV photons. 

Third, the flux density multiplied by the square of 
the distance to the source from an observer 
and the momentum rate in the outflow (dP/dt) of 
VLA la and VLA lb fit well with a known empirical 
relation (Anglada et al. 1992) and theoretical under- 
standing (Curiel et al. 1987) of the jet-induced model. 
Figure 2 shows the relation between S,^D^ and dP/dt for 
16 embedded objects, where we added VLA la/b with 
S^D'^ = 0.182 X 0.452 mJykpc^ (this work) and dP/dt = 
3 X 10"^ Mq yr-i kms'^ (Aso et al. 2000). Here, the lim- 
ited spatial resolution of the HCO+ and CO observations 
in Aso et al. (2000) did not resolve which of the three 
1.3mm clumps (MMS2-MMS4; Chini et al. 1997) is re- 
sponsible for the molecular outflow. However, VLA la 
and VLA lb in MMS2 are the only 3.6 cm sources that 
are associated with MMS 2-MMS 4 (this work; Reipurth 
et al. 1999). We therefore safely assumed that the dP/dt 
value determined for this molecular outflow represents the 
sum of the momentum rate from VLA la and VLA lb. 
Similarly, we summed the flux density of VLA la and 
VLA lb for the Si, value. 

Fourth, in case of VLA la, the source is elongated (the 
arrow in figure 1) in the direction of the global outflow 
seen in the H2 w = 1 -0 S(l)-band (Tsujimoto et al. 2002b; 
the main stream emanating from VLA la/lb in figure 3), 
which is characteristic of free-free centimeter emission in- 
duced by jets (Anglada 1996). 

3.2. The Origin of the Hard X-ray Emission 

The poor photon statistics of the hard X-ray source 
TKH 8 do not allow us to examine whether this is an ex- 
tended source. However, with the assumption that it is a 
point-like source, the X-ray emission is significantly offset 
from the class I protostars (IRS 3 and IRS 5; Tsujimoto 
et al. 2002b). Our high-resolution radio map enables us 
to confirm that the displacement is in the direction of 
the jet and outflow. We discuss the origin of this X-ray 
emission together with the centimeter emission based on 
a jet-induced scenario. 

Figure 4 shows a schematic view, where the hard X- 
rays are emitted from the post shock (PS) region, while 
the centimeter emission is from the recombination zone 
(RZ) behind the shock front. The RZ is maintained by the 
continuous ionization by UV photons from the PS region. 

In PS, the temperature (Tps) smd the density (npg) are 
expressed as 



Tps = 1.5 X 10'" 



100 kms" 




-2.5 -2 -1.5 -1 -0.5 
log {SJ)^ mJy kpc") 



0.5 



Fig. 2. Relation between Si/ and dP/dt. Open squares are 
from Anglada et al. (1992), who derived an empirical relation 
for these sources (dashed line). VLA la/b (filled square) is 
consistent with this relation. 
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Fig. 3. The H2 intensity map (gray scale) with the position 
of the 3.6 cm sources (pluses) and the hard X-ray intensity 
(contours). 
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(1) 
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nps = 4no [cm ^] , 



(2) 



where Vs is the velocity of the shock produced by the colli- 
sion of the jet into ambient matter with density no (Raga 
et al. 2002). When Vs is large enough, we can expect 
to have X-ray emission from PS. Assuming that light ele- 
ments are fully ionized in PS, the emission measure (EM) 
is given with the electron density (rips) and the volume 
(Vps) of PS as 

EM ^nlsVps[cm-% (3) 

We fitted the X-ray spectrum of TKH 8 with an optically- 
thin thermal plasma model with the amount of interstel- 
lar absorption, plasma temperature, and emission mea- 
sure as free parameters, obtaining an acceptable fit with 
Nh = 12.7 X 10^2 cm-2, fceTps = 3.05 keV, and EM = 
2.1 X 10^"^ cm^"^. Readers should note, however, that these 
parameters are not well constrained due to the poor X- 
ray photon statistics of this source, which can bring a 
large uncertainty in the following calculations. Using 
the best-fit values, we derive Vs = 1.5 x lO'^kms"^ and 
Uq = 5.8 X 10^ cm^^. Here, we assumed that PS is a cube 
at a distance of 450 pc with the length of 0."5 the pixel 
scale of ACTS). 

The radio observations give us another constraint on the 
values of Vs and uq as shown below. In RZ, the centimeter 
intensity (S^) is given by 



( = 1.42 X 10^ 
VmJyy 

lO^K 



^aOb 



arcsec^y v5GHz 
{1 -exp(-T^)}, 



(4) 



where 9a and Oh are the axis lengths of the extended cen- 
timeter source and Trz is the temperature in RZ (Curiel 
et al. 1993). The term 1 — exp(— r^) approaches in 
the optically-thin limit, which is expressed in terms of the 
shock parameters (Curiel et al. 1989) as 

1.68 



= 1.55 X 10" 



no 



5 GHz 



1cm 
-2.1 / r 



RZ 



104 K 



100 kms" 

-0.55 



(5) 



By substituting the observed values for VLA la 
(i/ =8.3 GHz, 5'^ = 0.128mJy, 9^ =0."40, and =0."09) as 
typical values and assuming that Trz = 10*^ K, we obtain 

,i£^)(iooiS^)'"-^-^^^«'- 

From the X-ray observations, we can independently de- 
rive that 

1.68 



Vlcm^sJ VlOOkms-i/ 



^ , , , , I =5.7x 1G\ (7) 
-Icm^^y VlOOkms-iy ^ ^ 

These values are in a good agreement with each other 
within a factor of a few, supporting the protostellar jet 
scenario for the origin of the X-ray and centimeter emis- 
sion. 
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Fig. 4. Schematic view of tlie centimeter and X-ray emission 
at the shock front. 



3.3. A Similar Example 

Similar X-ray emission was recently reported by Bally 
et al. (2003) from L1551 IRS 5, in which an X-ray source 
detected by Chandra is located at the base of the pro- 
tostellar jet. They argued three possibilities as the ori- 
gin of the X-rays, including X-rays from fast shocks. 
Applying our scenario to this source using the X-ray ob- 
servables given in Bally et al. (2003) and the centimeter 
observablcs in Cohen, Bieging, Schwartz (1982), we found 
that the discussion in Sect. 3.2 also holds for this source 
with Vs ~500kms~^. This velocity is reinforced by Pyo 
et al. (2002), who present NIR echelle spectroscopy re- 
sults on the [Fen] outflow from this source, and indicate 
the existence of a high velocity component reaching 400- 
500 km s-^ 

A recent series of works using [Fen] spectroscopy de- 
tected a highly-coUimated high velocity component ex- 
ceeding '^200kms^^ in protostar outflows, including DG 
Tau (250-300 km s"^ Pyo et al. 2003) and HL Tau 
(~250kms^^; T-S. Pyo, private communication). In op- 
tical jet observations, FS Tau B was measured both in 
radial velocity and proper motion, giving a jet velocity of 
-400kms-i (Eisloffcl & Mundt 1998). Future X-ray ob- 
servations of these systems with active jets are needed, to 
enrich the sample and examine the correlation between X- 
ray and centimeter obscrvables, in order to test our idea. 

3.4. Implications 

Protostellar jets at a speed of 500-1000 kms^^ may 
sound extreme. However, such a high velocity should not 
be surprising if we consider that even a Herbig-Haro ob- 
ject, which is far away from its powering source, is emit- 
ting X-rays induced by a ^200kms~^ shock (Pravdo et 
al. 2001), and also that solar coronal mass ejections are 
frequently observed to be propagating at >1000kms~^ 
(e.g., Gallagher et al. 2003). The lack of detection of 500- 
1000 km s~^ jets so far does not prove their non-existence 
but simply means we have had no tools to measure them. 

If this picture is established, we will obtain a strong 
method to derive the speed of highest velocity component 
of protostellar jets from X-ray imaging-spectroscopy using 
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equation (1), with no dependence on the inclination angle 
unlike proper motion or Doppler shift measurements. We 
consider more observations of similar sources will boost 
our understanding of this phenomenon, and eventually, of 
forming stars. 

4. Summary 

• We have made a high-resolution 3.6 cm imaging ob- 
servation of a hard X-ray emission source (TKH8) 
associated with a protostellar clump (MMS 2) in 
OMC-3 with the VLA. Two VLA sources were de- 
tected, both of which have the NIR counterparts. 

• The centimeter emission is concluded to be free-free 
emission produced by a protostellar jet based on 
(1) the association with class I protostars, (2) the 
low-mass nature of these protostars, (3) the relation 
between the centimeter flux density and the momen- 
tum rate of the outflow, and (4) the elongated struc- 
ture aligned with the global outflow. 

• The hard X-ray source is located ^ 1-2" offset in 
the direction of the protostellar jet and outflow from 
the NIR sources. The origin of the hard X-ray emis- 
sion was discussed in the framework of jet-induced 
shocks. The X-ray and the centimeter observations 
give independent constraints on the shock param- 
eters, which are found to be consistent with each 
other. 
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